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ABSTRACT 

We present the discovery of six millisecond pulsars (MSPs) in the High Time Reso- 
lution Universe (HTRU) survey for pulsars and fast transients carried out with the 
Parkes radio telescope. All six are in binary systems with approximately circular or- 
bits and are likely to have white dwarf companions. PSR J10T 7-7156 has a high flux 
density and a narrow pulse width, making it ideal for precision timing experiments. 
PSRs J1446-4701 and J1125-5825 are coincident with gamma-ray sources, and fold- 
ing the high-energy photons with the radio timing ephemeris shows evidence of pulsed 
gamma-ray emission. PSR J1502-6752 has a spin period of 26.7 ms, and its low period 
derivative implies that it is a recycled pulsar. The orbital parameters indicate it has 
a very low mass function, and therefore a companion mass much lower than usually 
expected for such a mildly recycled pulsar. 

In addition we present polarisation profiles for all 12 MSPs discovered in the 
HTRU survey to date. Similar to previous observations of MSPs, we find that many 
have large widths and a wide range of linear and circular polarisation fractions. Their 
polarisation profiles can be highly complex, and although the observed position angles 
often do not obey the rotating vector model, we present several examples of those that 
do. We speculate that the emission heights of MSPs are a substantial fraction of the 
light cylinder radius in order to explain broad emission profiles, which then naturally 
leads to a large number of cases where emission from both poles is observed. 

Key words: pulsars: general — pulsars: individual: (PSR J1017-7156, PSR J1337- 
6423, PSR J1446-4701, PSR J1502-6752, PSR J1543-5149, PSR J1622-6617) 



1 INTRODUCTION 

Radio pulsars are typically classified in terms of the observed 
parameters of pulse period and its derivative, placing them 
on the 'P — P diagram'. They are generally thought to be 
born with spin periods of a few tens of milliseconds and 
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magnetic field strengths of the order 10 G. As the pul- 
sar ages, the rotational energy is radiated away and the 
rotation rate slows, and therefore the pulsar traces a line 
across the P — P diagram, until it reaches the 'death line' 
where radio emission ceases. A distinct sub-group of radio 
pulsars, the "millisecond pulsars" (MSPs) typically have 
spin periods less than ~ 30 ms, often just a few millisec- 
onds. These periods are shorter than the youngest pulsars, 
however their small period derivatives indicate spin-down 
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ages of > 10 9 years and they are most often found in bi- 
naries with white dwarf (WD) or neutron star (NS) com- 
panions. The evolution of M SPs must therefore follow a dif- 
ferent path to the majority of pulsars. The accepted the- 
ory states that MSPs are initially formed in the same way 
as other pulsars, however after most are thought to have 
spun down to longer periods, they are then 'recycled' via 
accretion of mass and angu lar momentum from a compan- 
ion star (|Alpar et al.lll982 ). This allows the pulsar to at- 
tain very short spin periods (although this is not always the 
case), and the accretion process is thought to also reduce 
the dipole magnetic field str ength and therefore reduce the 
observed period derivative jBisnovatvi-Kogan fc Komberd 
ll974l ; lTaam fc van den Heuvelll 19861 ; ICumming et aLlbOOlT 
The majority of MSPs have spin periods less than 10 ms, 
and are either solitary or in circular orbits with a WD, or 
a very low mass companion (< O.O3M0). Longer spin pe- 
riod MSPs tend to have larger mass companions. Within 
10 < P < 20 ms the companions are always WDs, however 
above 20 ms there are a bro ad range of systems including 
WD. NS (iBurgay et aljliobl ) and even Be star companions 
(| Johnston et all 1 19921) At the time of writing, the pulsar 
catalogu^B lists 220 pulsars with a spin perioqf] less than 
30 ms, of which only 95 are located outside of globular clus- 
ters. Although the majority of the observed MSPs are lo- 
cated in globular clusters, the fact that many are formed 
via exchange interactions, and the effects of the cluster grav- 
itational potential means that their evolutionary paths are 
different. In this paper we only consider Galactic-field MSPs. 
Only 2% of pulsars with periods greater than 30 ms are in 
binaries, however there are 63 binary systems in the 95 non- 
cluster MSPs. 

The pulse profiles of the majority of MSPs share 
many characteristics with those of the classically 'young' 
pulsars , i.e. those with P < 100 ms and P > 
10" 15 (frhorsett fc Stinebrinelll99Ch . The profiles of MSPs 
are typically broa d with polarisation f raction between ~ 
1% and ~ 50% dXilouris et al l 1 19981 ; IStairs et al l 1 19991 ; 



lOrd et alj|2004 lYan et all 120111 ). One way to understand 



the polarisatio n of pulsars is through the 'r otating vector 
model' (RVM; iRadhakrishnan fc Cookdfl969l ) in which the 
observed polarisation angle is defined by the direction of a 
dipolar magnetic field at the point of emission. The model 
is defined in terms of the geometric parameters a, the an- 
gle between the rotation axis and the magnetic axis, and 
/?, the angle between the magnetic axis and line of sight 
at closest approach. For some pulsars, particularly those 
with components separated by ~ 180° this model has been 
used to me asure these geometric parameters with high pre - 
cision (e.g. iKramer fc Johnston! [20081 ; iKeith et af1l2010bl ). 
Although many MSPs have emission over a very wide phase 
window, sometimes greater than 180°, the position angles 
(PAs) often have large deviations from the best fit to the 
RVM, or t he RVM simply does not model t he observations 
well at all |Ord et alj|20o3 : lYan et al.ll201ll ). 

As well as the open questions intrinsic to the for- 



mation, evolution and emission of MSPs, they can also 
be valuable tools for other physics and astrophysics ex- 
periments. Recently many radio MSP s have been discov - 
ered to have gamma-ray pulsations |Abdo et al.l l2009bl ). 
and indeed searches of Fermi Large Area Telescope (LAT) 
gamma-ray sources have led to the dis covery of MSPs with 
both radio and gamma-ray emission (R ansom et al.l l201ll ; 
ICognard et~al]|201ll ; [keith et al.ll201ll ). This opens new pos- 
sibilities for studying the broad-band emission phenomenon 
of pulsar and comparison between young pulsars and MSPs. 
Although gamma-ray pulsations have been de tected from 
some young pulsars without rad io emission (|Abdo et al.l 
l2009al ; ISaz Parkinson et all |2010| ), the gamma rays from 
most MSPs can only be binned into pulse profiles with the 
aid of a timing ephemeris obtained through radio timing. 
The high precision time of arrival measurements that can 
be obtained from fast spinning MSPs, and their typically 
stable rotation rates mean that they can be used as precise 
clocks for tests of gravity and as broadband pulsed signals 
for probing the interstellar medium. Of particular interest 
is the use of a n array of MSPs to co nstrain the solar sys- 
tem ephemeris (|Champion et al]|2010l ) and for the detection 
of low-frequency gravitat ional waves jYardlev et alj l201ll ; 
Ivan Haasteren et al.|[201ll ). 

Research into the origin of the MSPs, the physics of 
their emission mechanism, and their use as tools for other as- 
trophysical experiments are all facilitated by expanding the 
known sample. To achieve this, the High Time Resolution 
Universe (HTRU) survey for pulsars and fast transients be- 
gan at the Parkes radio telescope in 2008. The high time and 
frequency resolution of the survey makes it mu ch more sen- 
sitive to MSPs than previous efforts at Parkes (|Keith et al.l 
l2010al ). I ndeed, we have already announced t he discovery of 
6 MSPs jBates et al.ll201ll ; lBailes et al.ll201ll ). In this paper 
we present the discovery and timing solutions for a further 
6 MSPs, and discuss their origins (Section [2J. We also anal- 
yse Fermi LAT data associated with PSRs J1125-5825 and 
J1406-4701, both of which are coincident with LAT sources 
(Section [3]). In Section 0] we present high quality polarised 
profiles for all 12 MSPs discovered by the HTRU survey to 
date and discuss the implications in Section [5] 



2 DISCOVERY OF SIX MSPS 

Here we report on the discovery of six MSPs in the HTRU 
survey, namely PSRs J1017-7156, J1337-6423, J1446-4701, 
J1502-6752, J1544-5149 and J1622-6617. All six MSPs are 
in almost circular orbits, consistent with the standard model 
in which the spin-up of the pulsar is associated with Roche 
Lobe overflow from a nearby companion that circularised the 
orbit. The companions are either white dwarfs or very low 
mass companions, such as t hose in orbit around PSR J2051- 
0827 (|Stappers et all 1 19961 ). We note that all 12 MSPs dis- 
covered by the HTRU survey to date are part of binary sys- 
tems, though there does not seem to be any obvious reason 
for the HTRU survey to preferentially select binaries. 



1 http: / /www. atnf.csiro.au/research/pulsar/psrcat/ 

2 Even though recycling process can produce pulsars with spin 
periods longer than 30 ms, we choose this limit to easily distin- 
guish MSPs from young pulsars. 



2.1 Observations and timing 

Each of these six pulsars was discovered in 540 s integrations 
as part of the mid-latitude portion of the HTRU survey be- 
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ing ca rried out at the Parkes radio telescope fsee lKeith et al.l 
l2010al for details). Follow up observations were carried out at 
Parkes with the centre beam of the 'HI-Multibeam' receiver 
and the Berkley-Parkes-Swinburne Recorder (BPSR), or 
once the orbital parameters were identified, the 3rd Parkes 
Digital Filterbank system (PDFB3). The BPSR observa- 
tions had an effective bandwidth of 340 MHz centred at 
1352 MHz and the PDFB3 observations had a bandwidth of 
256 MHz centred at 1369 MHz. PDFB3 observations were 
calibrated for differential gain between the linear feeds by 
means of observation of a pulsed noise diode. Each obser- 
vation was averaged in time, frequency and converted to 
Stokes total intensity, then matched with an analytic ref- 
erence profile to produce a single time-of-arrival measure- 
ment for each observation. A model of the rotational, astro- 
metric and orbital parameters was then fit to these mea- 
surements for each pulsar, in accordance with the stan- 
dard pulsar timing procedure, using the TEMPO 2 software 
(|Hobbs et alj|2006n . Due to the almost circular nature of the 
orbits we have used the 'ELL1' binary model as described in 
lLange et al.l (|200lT ). The best fit timing parameters and ad- 
ditional derived and observed properties of these six pulsars 
are presented in Tables [T] and [2] Flux densities presented 
are averaged over all 1369 MHz observations and are cali- 
brated by comparison to the continuum radio source 3C 218. 
We do not quote a standard error in this value since all of 
these pulsars exhibit significant flux density variations due 
to interstellar scintillation. Also shown in the ta ble is the oft 
used gamma-ray detectability measure E 1/2 /d 2 (|Abdo et al.1 
l201dh . MSPs detected by the Fer mi satellite typically have 
E 1/2 /d 2 > 10 10 erg 1/2 pc- 2 s- 1/2 (|Abdo et al.ll2009bf ). The 
measured frequency derivative can be biased du e to motion 
transverse to the line of sight (|Shklovskiill 1970l ) ■ No proper 
motion has yet been measured in these syste ms, however we 
can a ssume a typical velocity of 100 km s -1 (jToscano et al.l 
Il999l ) and the DM derived distance to compute a likely 
fraction of the frequency derivative which could be due to 
this effect. This is largest for PSR J1017-7156, with this 
value ~ 0.3 u, however the other five pulsars have a likely 
Shklovskii contribution less than 0.2 v. 



2.2 Two low mass binary pulsars: 

PSR J1017-7156 and PSR J1543-5149 

As seen in Tables [T] and H PSR J1017-7156 and J1543-5149 
have similar properties and fall into the category of 'low 
mass binary pulsars' ( LMBPs; IPhinnev fc Kulkarnil 1 19941 ; 
lEdwards fc Bailesll200ll ). Other than the low mass function 
(companion mass 0.1 < M c < 0.5 Mq), these pulsars are 
typified by an orbital period of a few days and a rapid ro- 
tation rate. It should be noted however that there are a 
number of so urces which fall into this category with longer 
spin periods (jEdwards fc Bailesl 120011 ). LMBPs are by far 
the most common class of binary MSPs, with some 30 (of 
57) such pulsars known. 

PSR J1017-7156 has a spin period of 2.34 ms and a DM 
of 94.2 cm -3 pc and was discovered in an observation centred 
~ 0.07° from the nominal pulsar position, carried out on 
2010-01-02 and detected with a signal-to-noise ratio (S/N) 
of 54. The candidate was initially id entified by a ne u ral ne t 
system similar to that described in lEatough et all (|2010h . 
Follow up observations quickly identified the pulsar to be 



in a 6.5 day circular orbit around a low mass companion. 
The properties of the binary are similar to those seen in 
other LMBPs. The high S/N, short pulse period and narrow 
profile suggest that this pulsar is a good candidate for high 
precision timing. Indeed, we can see that that this pulsar 
gives the best RMS residual of any HTRU discovery to date 
and it is alrea dy included as a ta rget for the Parkes Pulsar 
Timing Array jHobbs et al.ll2009l ). 

PSR J1017-7156 has a mean flux density of 0.89 mjy 
(the highest of any HTRU MSP), however even though its 
position was covered by previous surveys the relatively large 
DM ensured that it was not discernible from radio frequency 
interference due to excessive dispersion measure smearing. 
To demonstrate this , folding the closest pointing in the 
(|Edwards et al.l [20011 ) survey shows a signal with a S/N of 
11 and a pulse width of 1.3 ms which is broader than half 
the pulse period. The superior frequency resolution of the 
BPSR back-end allowed for an easy detection with a pulse 
width of 0.16 ms width and a S/N of 54. This clearly demon- 
strates the increased sensitivity of the BPSR back-end for 
pulsars with a high DM/P ratio. 

PSR J1543-5149 has a pulse period of 2.06 ms, orbital 
period of 8 days and was first detected with a S/N of 16 in 
an observation carried out on 2010-08-24, offset 0.08° from 
the nominal pulsar position. The minimum companion mass 
is 0.22 Mq , and so this pulsar is likely to be another member 
of the LMBP group. 



2.3 An intermediate mass binary pulsar: 
PSR J1337-6423 

Intermediate mass binary pulsars (IMBPs) have heavier 
companion masses than the LMBPs and typically have spin 
periods of order 10 ms. PSR J1337-6423 is a 9.4 ms pulsar 
with a DM of 260 cm -3 pc, the second largest DM of any 
MSP known to date. It is located 0.04° from the centre of 
the survey beam in which it was discovered with a S/N of 
11. The pulsar is in a 4.8 day orbit with an eccentricity of 
< 1 x 10~ 4 . These properties, and the minimum compan- 
ion mass of 0.74 M suggest that PSR J1337-6423 is likely 
to have a heavy WD companion and fall into the class of 
IMBPs. 



2.4 A very low mass binary pulsar: 
PSR J1446-4701 

PSR J1446-4701 is a 2.2 ms pulsar in a compact binary 
with a 6.6 hour orbital period. The pulsar was discov- 
ered in a survey beam centred 0.1° from the nominal pul- 
sar position on 2009-08-20 with a S/N of 12.4. The min- 
imum companion mass is just 0.019 Mq. With the recent 
discovery of PSRs J0610-2100 J1731 - 1845, J2214+3000 



and J2241-5236 jBurgav et al l 
iRansom et~aH |201ll ; iKeith et al.l 



20061; iBates et al.l 1201 ll ; 



201lh . all of which have 



very light compa nions, it is clear that the so-called "Black 
Widow" systems jKhizniak et al.lll988l ; IPhinnev et al.lll988l ; 
Ivan den Heuvel fc van Paradiisl Il988l ) are members of a 
larger class of v ery low mass binary pulsars (VLMBPs; 
iFreire et al.ll2003l ). These systems all show spin periods less 
than ~ 4 ms, orbital periods less than ~ 10 hours and a large 
E. Although many of these systems show eclipses for some 
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Table 1. Observed parameters from timing of three previously unknown MSPs. Values in parenthesis arc the nominal l-cr TEMP02 
uncertainties in the last digits. 



Parameter 



J1017-7156 



J1337-6423 



J1446-4701 



Right ascension, a (J2000) 10:17:51.32832(5) 13:37:31.928(18) 

Declination, 8 (J2000) -71:56:41.64003(12) -64:23:04.88(3) 

I (°) 291.56 307.89 

b (°) -12.55 -1.96 

Pulse frequency, v (s" 1 ) 427.62190510627(9) 106.11873502(3) 

Frequency derivative, v (s~ 2 ) -4.75(6) X 10" 16 -2.2(8) X 10" 15 

Epoch of model (MJD) 55329.1 55234.7 

Dispersion measure, DM (cm" 3 pc) 94.2256 260.32 

Binary model ELL1 ELL1 

Orbital period, P b (d) 6.5118988121(19) 4.78533407(9) 

Projected semi-major axis, asini (lt-s) . . . 4.83004527(11) 13.086499(11) 

Epoch of ascending node, T asc (MJD).... 55329.10065316(5) 55234.770356(6) 

ecosu -7.174(5)xl0" 5 1.82(13) x 10~ 5 

esinw 1.2268(5)xl0" 4 8.4(14)xl0~ 6 

Inferred eccentricity, e 1.4212(4) x 10~ 4 1.97(13) X 10~ 5 

Minimum companion mass, m Cjm i n (Mq) 0.19 0.78 

Fit time span (MJD) 55343.2—55681.2 55460.0—55694.3 

RMS of residuals (/its) 0.8 37.5 

Reduced x 2 1.3 1.7 

Mean flux density, S1400 (mjy) 0.89 0.32 

Pulse width at 50% of peak, W 50 (°) 10 23 

Pulse width at 10% of peak, W 10 (°) 20 52 

Spin down energy loss rate, E (ergs -1 ) . . 8.0 X 10 33 9.2 X 10 33 

Characteristic age, t c (years) 1.4 X 10 10 7.6 X 10 8 

Dipole magnetic field strength, B sur f (G) 7.8 x 10 7 1.4 x 10 9 

DM derived distance, d (kpc)* 3.0 5.1 

E 1 / 2 /d 2 {xl0 10 erg 1 /2 pc -2 s -i/2 ) 10 5 



14:46:35.71432(13) 

-47:01:26.7616(15) 

322.50 

11.43 

455.64401644177(16) 
-2.10(4)xl0 -15 
55647.8 
55.8340(11) 

ELL1 

0.2776660759(19) 

0.0640128(7) 

55647.8044387(7) 

1.3(18)xl0" 5 

l.l(21)xl0~ 5 

< 6 x 10~ 5 

0.019 

55358.6—55734.5 

2.5 

0.8 

0.37 

13 

48 

3.8 x 10 34 

3.4 x 10 9 

1.5 x 10 8 
1.5 

9.0 



of the orbit, as with PSR J0610-2100 and J2241-5236, PSR 
J1446-4701 does not show any sign of eclipse at 1369 MHz. 



2.5 Two intermediate-spin-period, low-mass 
binary pulsars: PSR J1502-6752 and 
PSR J1622-6617 

Although LMBPs typically have spin periods of a few ms, 
there are a handful of LMBPs that have spin periods 
greater than 10 ms. Examples are PSRs J1745-0952 and 
J1841 +0130, which have spin periods of 19 and 29 m s respec- 
tively ([Edwards fc Baile jEooil : lLorimer et al.ll2006t ). but are 
otherwise similar to the LMBPs. 

PSR J1622-6617, with pulse period 23.6 ms, perhaps 
fits into the category of intermediate period LMBPs. It was 
first observed on 2010-01-03, in a survey beam 0.07° from 
the nominal pulsar position and with a S/N of 12. We note 
that the observed spin and orbital parameters are similar 
to the known interm ediate spin period LMBP J1841+0130 
lLorimer et al.ll2006j ). 

A more puzzling situation presents itself with the dis- 
covery of PSR J1502-6752. Discovered with a S/N of 15.8 
in an observation taken on 2010-06-11 (offset 0.08° from the 
nominal pulsar position) its spin period is 26.7 ms and the 
mass function implies a minimum companion mass of just 
0.02 M , well below any known LMBP. Indeed, PSR J1502- 
6752 and J1622-6617 have the smallest mass functions of 



any of the binary 'intermediate period' MSPs (selected by 
P > 10 ms and P < 10 -18 ). The true companion mass does, 
of course, depend on the inclination angle, and it is possible 
that we are observing these LMBPs close to face on. For 
the companion mass of PSR J 1502-6752 to be greater than 
0.1 Mq, the inclination angle must be less than 13°, which 
has a probability of 0.025 of being drawn from a uniform dis- 
tribution of three-dimensiona l orientation s . It is also worth 
noting that, as mentioned bv lFreire et al.l (|2003l ). there is a 
distinct gap in the mass functions of the lowest mass binary 
MSPs between the VLMBPs, all of which have minimum 
masses less than 0.03 Mq and the LMBPs of which all have 
masses greater than 0.09 Mq. This suggests that there is a 
fundamental difference in the formation mechanism of the 
two systems. The gap suggests that the LMBPs with the 
lowest mass functions must be somewhat inclined to the 
line of sight, and therefore we feel that the anomalously low 
mass function of PSR J1502-6752 is unlikely to be due to 
an inclination effect. 

To better understand the unique nature of J1502-6752, 
we can consider where it lies amongst the known population 
of MSPs. By plotting the binary period and minimum com- 
panion mass of the known binary MSPs, as done in Figure 
[1] we can identify at least three main populations of binary 
MSPs. The VLMBPs appear in the lower left corner of the 
diagram, with masses lower than 0.03 Mq, and binary pe- 
riods which are typically shorter than a day. The LMBPs 
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Table 2. Observed parameters from timing of a further three previously unknown MSPs. Values in parenthesis are the nominal l-cr 
TEMP02 uncertainties in the last digits. 



Parameter 



J1502-6752 



J1543-5149 



J1622-6617 



Right ascension, a (J2000) 15:02:18.610(4) 

Declination, 8 (J2000) -67:52:16.78(2) 

I (°) 314.80 

b (°) -8.07 

Pulse frequency, v (s -1 ) 37.3909719910(4) 

Frequency derivative, v (s -2 ) — 4.0(3) X 10~ 16 

Epoch of model (MJD) 55421.2 

Dispersion measure, DM (cm - 3 pc) 151.75 

Binary model ELL1 

Orbital period, P b (d) 2.4844570(5) 

Projected semi-major axis, asini (lt-s) . . . 0.31756(3) 

Epoch of ascending node, T asc (MJD) .... 55421.21202(4) 

ecosu -5.3(125)xl0~ 5 

esinw -3.9(144) xl0~ 5 

Inferred eccentricity, e <2x 10~ 4 

Minimum companion mass, m CiIn i n (Mq) 0.022 

Fit time span (MJD) 55360.4—55757.5 

RMS of residuals (^ts) 67.9 

Reduced \ 2 0.7 

Mean flux density, S1400 (niJy) 0.68 

Pulse width at 50% of peak, W 50 (°) 40 

Pulse width at 10% of peak, W 10 (°) 

Spin down energy loss rate, E (ergs -1 ) . . 5.9 X 10 32 

Characteristic age, t c (years) 1.5 X 10 9 

Dipolc magnetic field strength, B aur f (G) 2.8 X 10 9 

DM derived distance, d (kpc)* 4.2 

E 1 / 2 /d 2 (xl0 10 erg 1 / 2 pc -2 s -1 /2) .2 



15:43:44.1498(10) 

-51:49:54.681(7) 

327.92 

2.48 

486.154232083(12) 

-3.8(2)xl0- 15 

55522 

50.92 

ELL1 

8.06077304(9) 

6.480281(5) 

54929.067833(7) 

2.02(11) xl0~ 5 

6.3(12)xl0- 6 

2.2(l)xl0~ 5 

0.22 

55432.3—55758.5 

13.5 

3.0 

0.70 
38 

34 



6.1 x 10' 
2.4 x 10 9 
1.7 x 10 s 
2.4 
4.2 



16:22:03.6669(6) 
-66:17:16.978(4) 
321.98 
-11.56 

42.33082901485(8) 

-1.14(4)xlO- 1(i 

55253.1 

87.94 

ELL1 

1.640635183(20) 
0.979380(6) 
55253.087284(3) 
4.3(117) xlO -6 
6.7(127) xl0~ 6 
< 2 x 10~ 5 
0.092 

55256.9—55733.5 

22.2 

1.3 

0.52 

50 

9.8 

1.9 x 10 32 
5.9 x 10 9 
1.2 x 10 9 
2.2 
0.3 
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Figure 1. Figure showing minimum companion mass and bi- 
nary period for all known binary pulsars with spin period below 
100 ms. Pulsars with discoveries presented in this paper are cir- 
cled. We have split the binary MSP population into three cate- 
gories as shown, the different symbols mark VLMBP (triangle), 
LMBP (square) and IMBPs (circle) as described in the text. 



contain the majority of the pulsars, and occupy a broad re- 
gion in the centre of the plot, showing a weak trend between 
companion mass and binary period. A third population, the 
IMBPs, have much more massive companions than LMBPs 
with similar binary periods. There is some overlap between 
the LMBPs and IMBPs at binary periods greater than a few 
days, however it is clear that there is a distinct population 
for short period binaries. We choose to split the LMBPs and 
IMBPs at by a minimum companion mass cut-off of 0.5 Mq. 
As described above, PSRs J1446-4701 and J1502-6752 fall 
into the VLMBP range, PSRs J1017-7156, J1543-5149 and 
J1622-6617 are LMBPs and PSR J1337-6423 is amongst 
those classified as IMBPs. 

Taking the same sample of pulsars, we now plot mini- 
mum companion mass against spin period, resulting in Fig- 
ure [5] In this plane the IMBPs and LMBPs are much harder 
to distinguish, but we retain the classification from Figure 
□ There are two clear outliers, the first is PSR J1903+0327, 
a known anomaly, which possesses a main sequence com- 
pani on, and is probably the descendant of a triple sys- 
tem (|Champion et al.ll20od iFreire et all 120111 ). The second 
clear outlier is PSR J1502-6752, which sits alone in a large 
space in the lower right corner of the figure. We must there- 
fore consider if PSR J1502-6752 formed through the same 
channel as the VLMBPs, or if it underwent some unique 
process. 



VLMBPs probably form when mass transfer begins be- 
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Figure 2. Figure showing spin period and minimum companion 
mass for all known binary pulsars with spin period below 100 ms. 
Pulsars with discoveries presented in this paper are circled. The 
different symbols mark VLMBP (triangle), LMBP (square) and 
IMBPs (circle) as described in the text. 
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Figure 3. Result of folding the LAT gamma-ray photons with 
the radio ephemeris of PSR J1446-4701 (histogram), overlaid 
with the 1369 MHz radio profile (curve). The profiles are aligned 



assuming a DM of 55.83 cm 
the radio flux density. 



1 pc, and an arbitrary scaling for 



fore a large degenerate core has formed as a result of the 
proximity of the neutron star to the companion. This is a 
stable process, and leads to a long period of spin-up and 
a low-mass companion that is often prone to ablation and 
given the right geometry, can eclipse the pulsar. How does 
PSR J1502— 6752 fit into this scenario? If we exclude the 
face-on orbit argument our only hint comes from the or- 
bital period of the system, which at 2.5 d is much greater 
than the other VLMBPs. Mass transfer in wider orbits hap- 
pens later in the donor star's life, leaving less time for mass 
transfer. What is unclear is why similar orbital period bina- 
ries have more evolved (and heavier) companions. So while 
the spin period of PSR J1502— 6752 is nicely explained by 
the larger orbital period, the mass of the companion is not. 
More speculative scenarios might involve the post-supernova 
collision of the neutron star with a still main sequenc e com- 
panion due to a favourab ly-oriented asymmetric kick (|Bailesl 
Il989l : iTauris et alJll999h and the subsequent disruption of 
the host star and mild spin-up of the pulsar. The formation 
mechanism of VLMBPs is not yet well understood, and the 
discovery of PSR J1502-6752 undoubtedly complicates the 
picture further, although the uniqueness of this system may 
indeed suggest a process that is rare and improbable. The 
discovery of further VLMBPs would be greatly beneficial in 
determining if PSR J1502— 6752 really is unique, or if the 
spin period distribution of the VLMBPs is in fact wider than 
currently understood. 



3 GAMMA-RAY ASSOCIATIONS 

Upon discovery of PSR J1446-4701 we determined that the 
gamma-ray source 2FGL J1446. 8-4701 lies only 2.7' from 
the nominal radio p osition, well with in the 95% confidence 
radius of r 95 = 7.2' (|Abdo et al.ll201ll ). From the radio tim- 
ing, we measure the pulsar's spin-down energy loss rate E = 

4 x 10 34 ergs -1 , and E 1/2 /d 2 = 9 x 10 10 erg 1/2 pc~ 2 s^ 1/2 , 
assuming a DM-derived distance of 1.5 kpc. This value is 
similar to that of ot her MSPs for whic h pulsed gamma rays 
have been detected (|Abdo et alj|2010h . 

Since the discovery of the HTRU MSP J 1125-5825 



ijBates et al.ll201lf ). the gamma-ray source 2FGL J 1125.0- 
5821 has been discovered, lying on ly 6' away, also w ithin 
the gamma-ray error radius, rqr,=9' (|Abdo et al1l201ll ). The 
energy loss rate of the pulsar is E = 7.9 x 10 34 and E 1/2 /d 2 = 
4 x 10 10 er g 1 ^ 2 pc~ 2 s - ^ 2 , as suming a DM-derived distance 
of 2.6 kpc ([Bates et al.ll201ll ). 

We searched for gamma-ray pulsations from both of 
these pulsars by phase binning Fermi LAT photons (arriving 
between MJD 54683 and 55774) using the radio ephemeris. 
We used "Pass-6 diffuse" class events (highest probability of 
being gamma-ray photons) and excluded events with zenith 
angles > 100° to reject atmospheric gamma rays from the 
Earth's limb. We used the bin-independent H-test statistic 
([de Jager fc Buschind[201ol ) to search over two parameters: 
the maximum angular separation from the pulsar position 
(0.1° < r < 2°), and the minimum photon energy cutoff 
(100 < iJcutoff < 1000 MeV), with a maximum energy fixed 
at 50 GeV. This truncates the point-spread function at low 
energies and decreases the number of background events. 
The best signal for PSR J1446-4701 (r = 0.8°, Setoff > 770 
MeV) is shown in Figure [3] with a post-trials significance of 
5.4-a, while the best signal for PSR J1125-5825 (r = 0.4°, 
-Ecutoft > 1000 MeV) is presented in Figure [4] with a post- 
trials significance of 4.9-ct. 

The alignment with the radio profile is accurate to 
about 0.01 in pulse phase as the DM is well constrained. 
The gamma-ray profile of PSR J1446-4701 is fairly broad 
and has its peak at a phase of 0.5 relative to the narrow radio 
peak. Although much lower in S/N, the radio and gamma- 
ray profile shape and alignment are somewhat similar to 
PSR J0437-4715 ([Abdo et al.ll2009bD . In PSR J 1125-5825 
the gamma-ray peak trails the radio peak by a phase of 0.6, 
which is large, but not uncommon for millisecond pulsars 
with complex radio profiles. The gamma-ray profile may be 
similar to PSR J2124-3358, which also shows a very com- 
plex radio profile, or perhaps similar to the double peaked 
profil es such as PSR J003 0+0451 but with an absent leading 
peak (|Abdo et al.ll2009bl ). 
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Table 3. Pulse period (P), DM and derived energy loss rate (E) obtained from radio timing of the 12 MSPs discovered in the HTRU 
survey to date. Note that E is not corrected for the Shklovskii effect, which may be large in some pulsars. Also provided is the rotation 
measure (RM). pulse width at 50% and 10% of the peak flux density (W50 and W10), as deriv ed from the 1369 MHz po larised pulse 
profile s described in this work. References for discovery and full parameters arc: [a] this work, [b] iBates et alj j2011T l and [c] iBailes et al.l 
l|201lh . 
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DM (cm~ 3 pc) 
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rgs -1 ) RM (radm -2 ) 
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2.34 


94.2 
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38 
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Figure 4. Result of folding the LAT gamma-ray photons with 
the radio ephemeris of PSR J1125-5825 (histogram), overlaid 
with the 1369 MHz radio profile (curve). The profiles are aligned 
assuming a DM of 124.79 cm -3 pc, and an arbitrary scaling for 
the radio flux density. 



4 POLARISATION PROFILES 

We now consider the polarisation propert ies of all 12 MSPs 
discovered in the HTRU survey to date (|Bates et al.ll20l"ll ; 
IBailes et al.ll20Tlh . the basic parameters of which are given 
in Table [3] To form high S/N profiles we summed all 
the observations used for timing of each pulsar with 64, 
256 and 1024 MHz wide bands centred at 732, 1369 and 
3100 MHz respectively. The signals were passed through a 
poly-phase filterbank and folded on-line using the Parkes 
Digital Filterbank System (PDFB3 and PDFB4). In the case 
of PSR J1017-7156 the ATNF-Parkes-Swinburne Recorder 
(APSR) was used to obtain coherently dedisperesd profiles 
at 732 and 1369 MHz. Observations at 732 and 3100 MHz 
were performed with the "10-50" receiver. Observations at 
1369 MHz were performed with the centre beam of the "HI- 
Multibeam" receiver. Before each observation a noise diode 
coupled to the receptors in the feed is observed to calibrate 
for differential gain and phase between the feeds. To cor- 
rect for cross-coupling of the receptors in the Multibeam 
receiver, we used a model of the Jones matrix for the re- 
ceiver computed by observation of the bright pulsar PSR 



J0437-4715 over the entire range of hour angles visible, using 
th e 'measurement equ ation modelling' technique described 
in Ivan Stratenl (|2004T ). The majority of timing is carried 
out at 1369 MHz, therefore profiles at the other frequencies 
are only included if the S/N ratio was high enough. Note 
that we follow the 'RVM sign convention' for the PA (see 
lEverett fc Weisberd l200ll ) . We measure the Faraday rota- 
tion observed towards each pulsar by fitting PA variations 
across the 256 MHz ba nd centred at 1 369 M Hz, using the 
algorithm described in iNoutsos et all (|200cf ). The best fit 
values of rotation measure (RM) are provided in Table 

We now present Figures [51 1161 showing profiles of the 
12 MSPs in total intensity, linear and circular polarisation. 
When multi-frequency data are present we include profiles at 
each of the frequencies, with an arbitrary phase alignment. 
Additionally we show the polarisation PAs, corrected with 
the RM from Table [3] and therefore absolute PAs as emit- 
ted at the pulsar, under the assumption that the RM used 
is correct. The figures are marked with a horizontal bar in- 
dicating the amount of uncorrected dispersion smearing due 
to the finite channel width of the DFBs. 

4.1 Description of profiles 

PSR J1017-7156 (Figure [5J. The profile of PSR 
J1017-7156 exhibits a single emission region with width at 
half the peak value of 10° at 1369 MHz, and a slight broad- 
ening at low frequencies. In general, th e frequency evolutio n 
of the profile is small, typical for MSPs (|Kramer et al.|[T999h . 
although there is some evolution noticeable on the leading 
edge between 732 and 1369 MHz. The linear polarisation 
fraction is about 30% at all three frequencies, however the 
circular polarisation goes from 20% at 732 MHz to 30% at 
1369 MHz to 45% at 3100 MHz. At 732 MHz the circular po- 
larisation is predominantly positive, however by 1369 MHz 
it is entirely negative in sign. Combined with the PAs (as- 
suming a rotation measure of —76 radm -2 ), this strongly 
suggests that the emission at 732 MHz is dominated by a 
polarisation mode orthogonal to that which is dominant at 
1369 and 3100 MHz. The polarisation PAs also show con- 
siderable variation with frequency, especially towards the 
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trailing edge of the profile. At 732 MHz the PAs show a 
downwards 'S-curve', however at f369 MHz the PAs exhibit 
a 'U' shape, and at 3f00 MHz, the PAs show an upwards 
slope. We discuss the atypical properties of the polarised 
emission observed from this pulsar in Section f5. II 

PSR 31125-5825 (Figure At 1369 MHz, the pro- 
file of PSR Jl 125-5825 exhibits three distinct components. 
Though the S/N is much lower, the profile shape is simi- 
lar at other frequencies, with some dispersion broadening 
apparent at low frequencies. At 1369 MHz the brightest of 
the three main components appears to be composed of at 
least three overlapping sub-components. At the higher fre- 
quency it appears that this component is starting to break 
up into three or four narrower peaks.. In total, emission ex- 
tends over ~ 200° of the profile. The brightest component 
exhibits varying linear and circular polarisation, with an or- 
thogonal mode jump aligned with the centre of the peak. 
There is also a change in the handedness of circular polari- 
sation associated with the jump. The other two components 
show no circular polarisation, one is 100% linearly polarised 
and the other is completely unpolarised. The spectral index 
of the components at phase —45° is flatter than that of the 
component at —80°. 

PSR J1337-6423 (Figure 0. The emission of PSR 
J1337-6423 extends over ~ 90° of pulse phase, with a half- 
width of 23°. This profile is ~ 5% polarised in both linear 
and circular, with the polarised intensity roughly following 
the shape of the total intensity. The PA shows a typical 
'S-shaped' swing. 

PSR J 1446-4701 (Figure 0). The profile of PSR 
J1446-4701 is fairly simple, showing a single peak with half 
width ~ 13° lead by a small shoulder. The profile has ~ 20% 
linear and circular polarisation, with the polarised intensity 
roughly following the shape of the total intensity. The po- 
larisation PAs are flat over the entire profile. The on the 
trailing edge of the peak (phases 0-20°), the ellipticity of 
the polarisation increases since the circular polarisation re- 
mains constant whilst the linear polarisation drops sharply. 

PSR J 1502-6752 (Figure The profile of PSR 
J1502-6752 covers nearly 300° of pulse phase. The main 
component has a polarisation fraction of ~ 40% linearly and 
~ 30% circular. The remainder of the profile does not have 
sufficient S/N to determine the polarisation fraction. The 
observed PAs appear to change linearly with pulse phase. 

PSR J1543-5149 (Figure [TOj). The profile of PSR 
J1543-5149 is composed of a single, broad, peak with small 
shoulders on the leading and trailing edges. There is less 
than 5% polarisation fraction for both linear and circular. 



45 




-0.2 - 



-0.4 - 

-0.6 h ',1',' 

-40 -30 -20 -10 10 20 30 40 

Pulse Phose (deg) 

Figure 5. Polarisation profiles of PSR J1017-7156 at (a) 732, (b) 
1369 and (c) 3100 MHz, showing phases within ±40° of the pulse 
peak. The solid black line shows total intensity, the thinner solid 
line shows linear polarisation and the dotted line shows circular 
polarisation. The inset figures show the profile over the full 360° of 
pulse phase. The data at 732 and 1369 MHz have been coherently 
dedispcrsed and so do not exhibit any DM smearing. 



PSR J 1622-66 17 (Figure El- The profile of PSR 
J1622-6617 is composed of a narrow peak with a broad 
shoulder on the trailing edge. The polarisation fraction is 
small throughout, however the linear and circular polarisa- 
tion track each other well, with a ~ 10% polarisation frac- 
tion in each. 
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Figure 9. Polarisation profile of PSR J1502-6752 at 1369 MHz. 
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Figure 10. Polarisation profile of PSR J1543-5149 at 1369 MHz. 
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Figure 11. Polarisation profile of PSR J1622-6617 at 1369 MHz. 



PSR J 1708-3506 (Figure [12]). The profile of PSR 
J1708-3506 appears to be formed of two main components, 
with the trailing component much broader than the lead- 
ing. The broad component appears to have a steeper spec- 
tral index, dominating the profile at 732 MHz and appearing 
only as a small shoulder at 3100 MHz, causing a relatively 
large profile evolution with frequency, as compared to most 
MSPs. The level of polarisation is less than 10% through- 
out, with the intensity of linear and circular following each 
other. An orthogonal mode jump is apparent ~ 5° after the 
peak of the profile. There appears to be a corresponding 
change in the handedness of circular polarisation associated 
with the jump, however this later reverses again without a 
corresponding change in PA. 

PSR J1719-U38 (Figure [13]). The profile of PSR 
J1719-1438 is relatively narrow, and shows a high shoul- 
der on the leading edge and a low shoulder on the trailing 
edge. The linear polarisation fraction is ~ 16% without de- 
tectable circular polarisation. The PA swing is steep with 
the steepest gradient leading the peak of the profile by 6° . 
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Figure 12. Polarisation profiles of PSR J1708-3506 at (a) 732, 
(b) 1369 and (c) 3100 MHz. 
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Figure 13. Polarisation profile of PSR J1719-1438 at 1369 MHz. 
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Figure 15. Polarisation profile of PSR J1801-3212 at 1369 MHz. 




of the circular polarisation may also reverse around, or just 
preceding this phase. 

PSR 3181 1-2404 (Figure HH}. At 1369 MHz we can 
discern at least four independent components with emis- 
sion spread over 250° of pulse phase. The emission has a 
mean linear polarisation fraction of 16% and a circular po- 
larisation fraction of 12%. Of the two brightest components, 
the leading component has a flatter spectral index and a 
higher polarisation fraction. The trailing component domi- 
nates the profile at 732 MHz, and the two components are 
equal at 3100 MHz. The profile appears to be almost com- 
pletely unpolarised at 732 MHz, even though the leading 
component has a considerable polarisation fraction at 1369 
and 3100 MHz. The inner components are also polarised and 
form a bridge of emission between the two main components. 



Figure 14. Polarisation profile of PSR J1731-1845 at 1369 MHz. 



PSR J 173 1-1 845 (Figure [14]). The profile of PSR 
J1731-1845 is composed of two main components, separated 
by 155°. The brightest component is polarised with a linear 
fraction of 50% and circular fraction of 20%. The second 
component is much less polarised, with only a hint of linear 
polarisation. 

PSR 31 80 1-321 2 (Figure [15]). The profile of PSR 
J1801-3212 consists of a bright leading component with a 
shoulder. The emission is elliptically polarised with a linear 
and circular polarisation fraction of ~ 15% each. The peak 
of the polarised intensity trails the peak of the total intensity 
by ~ 18°. The PA does not appear to follow the typical 'S- 
shaped' swing predicted by the RVM. The emission changes 
hand of circular polarisation at a phase corresponding to the 
peak of the total intensity profile. Just prior to this phase, 
the PA of the linear polarisation jumps by ~ 50°. The lim- 
ited S/N makes it hard to discern, but there may also be a 
true 90° jump in the PA between the polarisation associated 
with the main peak and that of the shoulder, occurring at a 
phase of ~ 40°. There may also be a hint that handedness 



5 DISCUSSION OF POLARISATION AND 
SPECULATION ON GEOMETRY 

We can use the RVM to model the observed polarisation PA 
>]/ in terms of geometrical parameters a and fi by 



tan(* + * ) 



sin a sua 



sin£ cos a — cos C, sin a cos{4> — (j>o) ' 



(1) 



where £ = a + f3, and ^0 and <f>o are the fiducial phase and 
PA, nominally corresponding to the magnetic axis crossing 
the line of sight. The relation between <j>o and the observed 
emission is however affected by aberration and retardation 
effects in the pulsar magnetosphere. If the emission height, 
r cm , is small compared to the light cylinder radius, -Rlc, 
then <f>o arrives later with respect to the corresponding total 
intensity emission, with a shift given by 

= 4r crn /R LC (2) 

(iBlaskiewicz et al.l Il99ll : iHibschman fc Aronsl l200ll ; iDvkd 
l200Sf ). 

The magnetospheres of MSPs are very compact, a pul- 
sar with a spin period of 5 ms has a light cylinder radius 
of only 240 km. Conventionally it is assumed that the radio 
emission occurs from a radius which is only a small fraction 
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Figure 16. Polarisation profiles of PSR J1811-2404 at (a) 732, 
(b) 1369 and (c) 3100 MHz. 



(few percent) of the light cylinder alt hough in young pulsars 
highe r emission heights are common (jjohnston fc Weisberd 
l200rJ ). In MSPs, low emission heights cause a problem be- 
cause an emission height of 10% of the light cylinder yields 
a pulse width of only ~60° unless the rotation and mag- 
netic axes are close to alignment. Observationally it has 
been clear for over a decad e that MSP widths are much 
large r than this on average (jXilouris et al.|[l998l ; lYan et all 
l201jjh . Indeed, the observed large widths preclude any emis- 
sion height less than the light cylinder in the conventional 
picture, unless the majority of MSPs are aligned rotators. 
Conversely, for a number of MSPs where there has been an 
attempt to constrain the geometry, it has been found that 
the inferred opening angle is much smaller than that ex- 
pected by extrapolation fro m the main population of radio 
pulsars (|Kramer et al.|[l998l ). 

In principle, the polarisation properties could help shed 
light on this conundrum, especially as large profile widths 
are a boon to RVM fitting. Unfortunately, the PA variations 
in MSPs can be extremely complex and are often far from 
that expected in the RVM. Furthermore, the lack of obvious 
inflexion point in the PA variations and the lack of identi- 
fiable core/cone components makes it difficult to estimate 
emission heights via the aberration/retardation methods of 
iBlaskiewicz et al] (|l99ll ). We note in passing however that 
MSPs appear to be no better or worse than normal pul- 
sars in this regard (lEverett fc Weisberd l20o"ll ) . Indeed the 
geometry of the wide profile pulsars B0950+08, B1822-09 
and B1929+10 to take but three examples remains unclear 
ijEverett fc Weisberd lioOll ). In spite of these inherent diffi- 
culties we have attempted RVM fits to several of the MSPs 
in our sample in an attempt to get a handle on geometry 
and emission locations. 

5.1 PSR J1017-7156 

Neither the drastic change of polarisation with frequency, 
nor the 'U' shaped PAs observed in this pulsar can be sim- 
ply explained within the framework of the RVM. We also 
observe a change in the handedness of circular polarisation, 
and a ~ 90° change in the absolute PA of the linear po- 
larisation between the profile at 732 MHz and 1369 MHz. 
We have considered the possibility that this effect is caused 
by an error in calibration, however multiple observations of 
this pulsar and the lack of similar effects in other pulsars 
largely rules this out. This is, therefore, a strong suggestion 
that the profile is composed of two competing elliptically 
polarised orthogon al modes, with different spectral indices 
(ISmits et al ] |2006h . This superposition of modes can also ex- 
plain the rapid sweeps in PA observed at the trailing edge 
of the pulse at 732 and 1369 MHz. The strong depolarisa- 
tion of the profile at the trailing edge suggests that there is 
a change in the dominant polarisation mode at this phase, 
and the PA should jump by 90°. Broadening of the pulse 
in the interstellar medium can smooth these discrete jumps 
somewhat, causing the appearance of a rapid swing in PA 
( Nou tsos et al. I l2009l ; lKarastergioull2009h . 

5.2 PSR J1125-5825 

The 'main pulse' of this pulsar extends over 100° of lon- 
gitude and has a well behaved swing of polarisation once 
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Figure 17. Polarisation profile of PSR J1125-5825 at 1369 MHz, 
showing the observed PAs overplayed with the model obtained 
from fitting the RVM. The dashed line is separated by 90° from 
the model PA to indicate orthogonal modes. 

the orthogonal jump near the centre of the profile is taken 
into account. The presence of a further polarised component 
some 130° away from the main pulse is an added complica- 
tion; it is not possible to fit the RVM to all the data. How- 
ever, it is possible to obtain a good RVM solution using only 
the main pulse data. In this case, a — 128°, /3 = —21° with 
the inflexion point coincident with the negative peak of the 
circular polarisation. The RVM correctly predicts the swing 
of the polarisation at the distant component but the values 
of PA are ~30° offset from the observed values. The data 
and model are shown in Figure [TT) 

We note also the similarity b etween this profile 
and that of the MSP J1012+5307 (jXilouris et alj 1 19981 : 
IStairs et al]|l999h and also wit h the normal pulsar B1055- 
52 |Weltevrede fc Wright]|2009h . In the former case the au- 
thors were unsure whether the pulsar is an aligned or an 
orthogonal rotator whereas the latter appears to be clearly 
orthogonal. For PSR J1125-5825 the ambiguity remains due 
to the problematical nature of the RVM fit. If the pulsar is 
an aligned rotator this could explain the large width and 
the complex RVM can arise if different components emit at 
different heights. On the other hand, if the pulsar is an or- 
thogonal rotator the likely implication is that the main and 
interpulse emission arise from significantly different heights 
or indeed from the outer magnetosphere, resulting in a non- 
standard RVM. With more time, and therefore greater S/N, 
the gamma-ray profile of PSR J1125-5825 may be able to 
provide further constraints on the geometry, and clues to 
the location of the radio emission regions. 

5.3 PSR 1502-6748 

This pulsar appears to be a very wide triple profile with the 
central peak dominating, and two weak outriders located 
some 80° away. The PA swing is smooth, the RVM fit is not 
very constraining on a and /3 but does locate the inflexion 
point some 40° distant from the main peak. This yields a 
height of 220 km, about 20% of the light cylinder radius, 
using Equation[2] Although this seems plausible it is difficult 



to reconcile this relatively low height with the large pulse 
width unless a is small. 

5.4 PSR J1708-3506 

The emission from this pulsar seems to conform to standard 
behaviour. It is a blended triple, with the central component 
having a steep spectral index and located equidistant from 
the flatter spectrum conal components. The PA swing is rel- 
atively smooth apart from an orthogonal jump. The RVM 
fit is reasonably well constrained with a between 60° and 
90° and /3 of 20°. The inflexion point occurs some 40° after 
the location of the central component. If we interpret the 
offset as due to aberration, this yields a height of ~40 km, 
which, assuming an emission region filled to the edge of the 
open field line region, would imply a profile width of ~80°. 
The observed profile is however much broader than this, 
placing doubt on the RVM-derived geometry, or otherwise 
suggesting a more complex model of the emission region is 
required. The large frequency-dependant evolution of the 
profile may suggest that the profile is composed of compo- 
nents that are originating in different parts of the magneto- 
sphere and therefore confusing our picture. 

5.5 PSR J1719-1438 

This pulsar exhibits an almost classical S-shaped swing of 
linear polarisation across the pulse. The profile is again ex- 
tremely wide, and this width coupled with the PA swing 
implies a small value of both a and /3. The inflexion point 
of the traverse is located some 6° prior to the pulse peak. 
The geometry of this star is thus rather uncertain and can- 
not shed much light as to the origins and ev olution of its 
planetary mass companion jBailes et alj|201ll ). 

5.6 PSR J1811-2404 

Polarised emission extends over nearly 300° of longitude in 
four distinct components. By inserting an orthogonal mode 
jump between the dominant component and the weaker 
components we can obtain a remarkably good and strongly 
constrained RVM fit with a=89.7, /3=21 and the location of 
the inflexion point 145° after the main pulse. The data and 
model are shown in Figure [TH] The pulsar therefore seems to 
be an orthogonal rotator where we see emission from both 
poles. The question remains as to which components belong 
to which pole, and which (if any) are centrally located. There 
are two possibilities. 

The first option is that the core is located at phase 
90° and —90°. The main pulse would then consist only of 
the core component with no conal outriders. The interpulse 
consists of the core plus two trailing conal components. The 
spectral index evolution does not fit this picture very well, 
and the extreme width of the interpulse is also difficult to 
reconcile. An alternative form of this picture would assign 
the two central components to a different emission height, or 
from an emission region located in the outer magnetosphere. 

In the second possibility, the core is located at phase 
135° and —45° . In this case, the main pulse consists only of 
a leading edge cone whereas the interpulse is a mostly sym- 
metric triple structure. Under this picture the offset between 
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Figure 18. Polarisation profile of PSR J1811-2404 at 1369 MHz, 
showing the observed PAs overlaid with the model obtained from 
fitting the RVM. The dashed line is separated by 90° from the 
model PA to indicate orthogonal modes. 



the inflexion point and the core location yields a height of 
54 km (c.f. the light cylinder radius of 124 km) and this 
height is consistent with a pulse width of some 120°, about 
that of the interpulse emission. 



5.7 Summary 

There now exists a substantial body of high quality polari- 
sation measurements of MSPs built up over the last decade. 
Key points to re-iterate from these observations are (a) 
many MSPs have lar ge widths especia lly when observed with 
high dynamic range (|Yan et al.ll201ll ), (b) as in normal pul- 
sars a wide range of linear and circular polarisation fractions 
are observed (c) the polarisation profiles can be highly com- 
plex and the PA swing often does not obey RVM, (d) a 
non-negligible fraction are amenable to RVM fitting. 

It seems difficult to reconcile the observations with a 
low emission height assuming an isotropic distribution of 
cos(a). However, it may be that a is biased towards small 
values in MSPs, linked to their formation process. This has 
to be investigated further but the presence of bridges of 
emission between the various components might support this 
idea. However, it is also not a given from the data that the 
entire phenomenological view of the radio emission has to 
be overturned in MSPs by invoki ng e.g. caustic em ission 
or extreme magnetosphere effects (|Dvks et alj|2010l ). Con- 
sider a range of allowable emission heights up to ~50% of 
the light cylinder radius. For isotropic values of cos(a) one 
would then expect minimum pulse widths of ~50° and high 
fraction (30%) of interpulses. This is more or less consis- 
tent with the data where it could be a r gued that at least 3 
of our sample and 7 of the lYan et aLI (|201ll ) sample show 
evidence of interpulse emission. In these interpulse pulsars 
there is no a priori reason why either the flux density of 
the emission from the two poles should be similar (e.g the 
extreme case of PSR J1603— 7202) nor that the emission 
heights should be similar which would lead to highly com- 
plex PA traverses. The bridge emission seen in e.g. PSRs 
1939+2134 and 1045-4509 is not evidently part of the in- 



terpulse picture but the wide beams mean that naturally 
there will be blending between the main and interpulse pro- 
files, so we do not see this as a major issue. For the pulsars 
which clearly do not show interpulse emission, some have 
relatively narrow profiles (PSRs 1017-7156, 1022+1001, 
1446-4701, 1600-3053, 1713+0747) whereas some are wide 
(PSRs 1543-5149, 1708-3506, 1732-5049) perhaps indicat- 
ing lower values of a. In summary, we believe that emission 
heights in MSPs are high, up to 50% of the light cylinder. 
This leads to a substantial fraction of MSPs with interpulses, 
further complicating the identification of components in the 
profile and the PA swing. 

Although some MSPs have polarisation consistent with 
the RVM, it is still not clear that the derived geometry can 
be used to determine a consistent model for the emission lo- 
cation, and it is possible that the accretion process itself has 
modified the magnetosphere to a point where a simp le dipo- 
lar model is no longer applicable (|Xilouris et al.|[T998h . Over- 
all, however, we find that the polarised emission of MSPs is 
not more complex than that observed in many slow pul- 
sars. There has been success in untangling these complex 
PA variations in slow pulsars through st udies of polarised 
inten s ity on a pu l se-by - pulse basis (e.g. Backer fc RankirJ 
1 19801 ; iGil fc Lvnel Il995l ; iKarastergiou et al.l |2002| ). Future 
high-sensitivity observations may well be able to repeat this 
success in MSPs, giving us greater confidence in geometric 
interpretations of the RVM and a fuller picture MSP emis- 
sion. High significance gamma-ray profiles of MSPs could 
also provide additional constraints on geometry, and test 
theories of radio emission from the outer magnetosphere. 



6 CONCLUSION 

The High Time Resolution Universe survey for pulsars and 
fast transients has discovered 12 MSPs to date, six of which 
are announced in this work. These six MSPs (indeed, all 12) 
are in binaries. Two, PSRs J1017-7156 and J1543-5149, with 
spin periods of ~ 2 ms, binary periods of about 1 week and 
companions with m c ,mm ~ 0.2 Mq are easily classified as 
LMBPs, the most common category of MSPs. PSR J1337- 
6423, with its larger companion mass (m Cjm i n = 0.8 Mq) 
and longer spin period (9.4 ms) better fits into the class 
of IMBPs, likely with a heavy WD companion. In con- 
trast, PSR J1446-4701 has a much lower mass companion 
(nic.min = 0.02), falling easily into the distinct category of 
VLMBPs, typified by the 'black widow' pulsars. The high 
rotational energy loss rate, E ~ 4 x 10 34 , and inferred dis- 
tance of only 1.5 kpc suggests that this is a good candidate 
for detection as a gamma-ray pulsar. Indeed, we find a tanta- 
lising indication of gamma-ray pulsations by using the radio 
ephemeris to fold gamma-ray photons detected by the Fermi 
LAT. 

The final two pulsars, PSRs J1502-6752 and J1622- 
6617, both have spin periods of ~ 25 ms, typically considered 
intermediate spin periods for recycled pulsars. PSR J1622- 
6617 can be classified as an intermediate period LMBP, sim- 
ilar to other known systems such as PSR J1841+0130. There 
are however no similar systems to PSR J 1502-6752 in the 
literature. The mass function implies m Cjm i n = 0.02, typical 
of the VLMBPs, however the formation of these systems is 
thought to spin the pulsar up to very short periods. Whilst 
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it is possible that the low mass function is an inclination 
effect, it is certainly striking that there are no MSP systems 
with 0.026 < Tti'C, min < 0.09, with a large number of MSPs 
having m c ,mm ~ 0.1 or m c ,min ~ 0.02. Therefore we con- 
clude that the companion of PSR J1502-6752 is likely to 
have formed through the same channel as the companions 
for the VLMBPs, and so that formation mechanism does not 
require, or does not guarantee, a short spin period MSP. 

We have also undertaken polarimetric observations of 
all 12 of the MSPs discovered by the HTRU survey to date. 
The calibrated profiles show the wide variety of profiles typ- 
ical of MSPs, with pulse widths ranging from 24° to 280° 
(measured at 10% of peak flux). The profiles generally con- 
sist of multiple components, often showing broad shoulders 
running off either the leading or trailing edge of a narrow 
central feature. PSRs J1125-5825, J1731-1855 and J1811- 
2404 show the dominant components separated by ~ 180°, 
and may either be classified as interpulse pulsars or ex- 
tremely wide cones. For a number of the MSPs, we find 
that the observed swing in PA can be fit by the RVM, how- 
ever there are clearly cases such as PSR J1017-7156 where 
effects such as smoothed orthogonal mode jumps make this 
impossible. We believe that emission heights in MSPs are a 
substantial fraction of the light cylinder, leading to a large 
fraction of MSPs showing emission from both poles. 
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